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Electrical  Transport  Properties 
of  Polycrystalline  Monolayer 
Molybdenum  Disulfide 

Sina  Najmaei,t§  Matin  Amani/,s  Matthew  L.  Chin/  Zheng  Liu,T"  ¥  A.  Glen  Birdwell/ 

Terrance  P.  O'Regan/  Pulickel  M.  Ajayan/  Madan  Dubey,**  and  Jun  Lour  * 

department  of  Materials  Science  and  NanoEngineering,  Rice  University,  Houston,  Texas  77005,  United  States  and  ‘Sensors  and  Electron  Devices  Directorate, 
US  Army  Research  Laboratory,  Adelphi,  Maryland  20723,  United  States,  dhese  authors  contributed  equally.  Present  address:  School  of  Materials  Science  and 
Engineering,  Nanyang  Technological  University,  639798,  Singapore.  “Present  address:  NOVITAS,  Nanoelectronics  Centre  of  Excellence,  School  of  Electrical  and 
Electronic  Engineering,  Nanyang  Technological  University,  639798,  Singapore. 


ABSTRACT  Semiconducting  MoS2  monolayers  have  shown  many  promising 
electrical  properties,  and  the  inevitable  polycrystallinity  in  synthetic,  large-area 
films  renders  understanding  the  effect  of  structural  defects,  such  as  grain 
boundaries  (GBs,  or  line-defects  in  two-dimensional  materials),  essential.  In  this 
work,  we  first  examine  the  role  of  GBs  in  the  electrical-transport  properties  of  MoS2 
monolayers  with  varying  line-defect  densities.  We  reveal  a  systematic  degradation 
of  electrical  characteristics  as  the  line-defect  density  increases.  The  two  common 
MoS2  GB  types  and  their  specific  roles  are  further  examined,  and  we  find  that  only 
tilt  GBs  have  a  considerable  effect  on  the  MoS2  electrical  properties.  By  examining  the  electronic  states  and  sources  of  disorder  using  temperature- 
dependent  transport  studies,  we  adopt  the  Anderson  model  for  disordered  systems  to  explain  the  observed  transport  behaviors  in  different  temperature 
regimes.  Our  results  elucidate  the  roles  played  by  GBs  in  different  scenarios  and  give  insights  into  their  underlying  scattering  mechanisms. 

KEYWORDS:  two  dimensional  materials  •  MoS2  •  polycrystallinity  •  grain  boundaries  •  electronic  transport 


The  existence  of  an  intrinsic  direct 
band  gap  and  the  fascinating  physics 
of  two-dimensional  (2D)  MoS2  have 
led  to  intensive  studies  of  this  material  in 
the  past  few  years.'  7  This  effort  has  been 
greatly  facilitated  by  the  development  of 
new  MoS2  synthetic  approaches.8  12  The 
introduction  of  a  large  number  and  diverse 
types  of  defects  into  the  crystal  struc¬ 
ture  of  this  material  is  a  natural  consequence 
of  chemical  growth."  14  Therefore,  under¬ 
standing  the  nature  of  these  disorders  and 
their  roles  in  the  physical  properties  of  MoS2 
is  essential  for  its  future  development.  A 
grain  boundary  (GB),  a  prominent  type  of 
defect  in  chemical  vapor  deposition  (CVD)- 
grown  MoS2,  is  a  one-dimensional  (ID)  line 
defect  in  the  2D  lattice  and  mediates  the 
transition  from  one  crystalline  orientation 
to  another.  Previous  research  on  the  role 
of  grain  boundaries  (GBs)  in  the  transport 
properties  of  2D  materials  shows  consider¬ 
able  diversity.14-18  The  theoretical  examina¬ 
tion  of  their  roles  in  graphene  predicts  two 


distinct  behaviors  based  on  different  grain 
boundary  structures.15  Li  et  at.  demonstrated 
the  collective  role  of  grain  boundaries  in  the 
transport  properties  of  graphene  by  control¬ 
ling  the  grain  size  from  6  to  20  y/m.17  These 
results  show  that  the  charge  carrier  mobility 
of  the  samples  decreases  roughly  by  2-fold 
as  the  grain  sizes  decrease  in  this  range. 
Additional  experiments  on  graphene  samples 
synthesized  through  different  synthetic  routes 
suggest  that  sample  preparation  may  play  an 
even  more  important  role  in  electron  transport 
behaviors  across  grain  boundaries.18  These 
results  show  that  growth  conditions  may  dom¬ 
inate  the  quality  of  grain  boundaries  formed. 
For  instance,  slower  graphene  growth  results  in 
larger  grain  sizes;  however,  the  layers  are  often 
discontinuous  across  grain  boundaries.  For  this 
reason,  electron  transport  across  grains  made 
from  this  growth  process  show  higher  degra¬ 
dation  in  transport  properties  of  graphene. 

The  polycrystalline  MoS2  films  show  simi¬ 
lar  complexity  in  GB  structures,  and  dif¬ 
ferent  forms  of  line-defects  have  been 
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observed.11-14  It  has  been  shown  that,  in  addition  to 
conventional  GBs,  bilayer  growth  with  imperfect  bond¬ 
ing  is  also  possible.11  Additionally,  since  the  kinetics  of 
nucleation  and  growth  in  CVD-grown  MoS2  is  diverse, 
examining  the  effects  of  GBs  in  large-area  polycrystal¬ 
line  samples  in  different  scenarios  is  essential  to  fully 
comprehend  their  role.  Recent  reports  on  the  electron 
transport  characteristics  of  individual  grain  boundaries 
in  monolayer  MoS2  by  examining  merged  single¬ 
crystalline  domains  have  demonstrated  some  interest¬ 
ing  results,12  but  detailed  measurements  on  more 
diverse  sets  of  grain  boundaries  and  understanding  of 
the  underlying  scattering  mechanisms  involved  in  these 
GBs  remain  elusive. 

Here,  we  explore  field-effect  transport  characteris¬ 
tics  of  CVD-grown  MoS2  with  a  focus  on  the  role  of  GBs. 
We  study  the  field-effect  characteristics  of  polycrystal¬ 
line  large-area  MoS2  films  to  elucidate  the  collective 
role  of  GBs  in  the  electron  transport.  To  gain  insight 
into  these  results,  we  further  examine  two  common 
MoS2  GB  types  and  determine  that  only  tilt  boundaries 
have  a  degrading  effect  on  the  MoS2  electrical  transport 
properties.  Through  temperature-dependent  studies  of 

a 


the  GBs,  we  explore  the  mechanisms  of  charge  carrier 
transport  and  localizations  and  provide  a  comprehen¬ 
sive  understanding  of  the  GB  related  transport  charac¬ 
teristics  of  monolayer  MoS2. 

RESULTS 

To  study  the  collective  effect  of  GBs  on  the  electrical 
performance  of  CVD-grown  MoS2,  one  must  control 
and  quantify  the  density  of  line  defects.  Li  et  al.  de¬ 
signed  such  experiments  based  on  controlling  the 
grain  sizes  in  the  films.17  A  systematic  control  of  grain 
sizes  in  MoS2  is  notfully  developed,  and  consequently, 
we  introduce  two  alternative  strategies  for  control  of 
I  ine-defect  densities:  strategy  1 ,  systematic  variation  of 
the  line-defect  densities  in  devices  with  varied  channel 
lengths;  strategy  2,  stochastic  variation  of  line-defect 
densities  in  devices  with  similar  channel  lengths.  We 
utilize  these  two  complementary  strategies  to  examine 
the  statistical  electrical  transport  behaviors  in  poly¬ 
crystalline  samples.  For  both  of  these  approaches,  one 
needs  to  develop  a  line-defect  density  quantification 
approach  and  substantiate  the  controlled  nature  of 
changes  in  the  line-defect  densities.  A  common  feature 

b 


g 


Channel  Length  (pm) 


Figure  1.  (a)  Large  area  films  suitable  for  study  of  collective  role  of  grain  boundaries.  The  arrows  highlight  instances 
of  overlay  growth,  commonly  seen  at  grain  boundaries,  (b)  Close  examination  of  secondary  layer  growth,  a  common  feature 
of  grain  boundaries,  using  dark  field  transmission  electron  microscopy,  (c,  d)  The  SEM  and  optical  image  of  the  device 
designed  for  the  channel  length  dependency  measurements,  (e)  Corresponding  Raman  map  of  the  same  location,  showing 
the  distributions  of  the  grain  boundaries  in  the  film.  <f)  We  estimate  the  line  defect  density  by  calculating  the  total  length  of 
the  bilayered  regions  (25.5/rm  in  this  example)  and  dividing  by  the  total  channel  area  (1 80  /rm2  in  this  example),  resulting  in  a 
line  defect  density  of  0.14/rm  ’  in  the  MoS2  ribbon,  (g)  The  channel  length  dependency  of  the  estimated  line  defect  density. 
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Figure  2.  (a)  Field  effect  characteristics  of  devices  with  similar  cannel  length  of  1 8/im  but  with  varied  estimated  line  defect 
densities  containing  a  variety  of  randomly  oriented  and  grain  boundary  types.  Inset  represents  the  calculated  mobilities 
for  these  devices  as  a  function  of  line  defect  density,  (b,  c)  Estimated  changes  in  the  mobility  and  on/off  ratio  of  devices  as 
a  function  of  channel  length. 


of  large-area  MoS2  films  as  confirmed  by  electron 
microscopy  is  the  growth  of  secondary  layers  along 
GBs  (Figure  1a,b).  Many  macroscopic  tools  such  as 
scanning  electron,  optical,  or  Raman  microcopy  can 
be  used  to  identify  these  features  (Figure  1c  e).'9 
Assuming  that  the  formation  of  secondary  layer  stripes 
at  the  grain  boundaries  is  uniform  across  the  films, 
one  can  use  this  easily  quantifiable  characteristic  as  an 
indirect  measure  to  acquire  a  rough  estimation  of  line- 
defect  density  in  the  films.  We  approximate  the  line- 
defect  density  by  measuring  the  total  length  of  these 
observable  line  defects,  normalized  by  the  total  area  of 
their  channels  (Figure  If). 

For  strategy  1,  we  fabricated  seven  devices  with 
varied  channel  lengths  and  characterized  their  line 
defect  densities.  In  each  device,  the  large-area  MoS2 
films  were  patterned  into  ribbons,  and  FET  transistors 
with  different  channel  lengths  were  prepared  using 
conventional  photolithography  and  metallization  pro¬ 
cesses.  We  observe  a  systematic  change  in  average  line 
defect  density  of  these  devices  as  a  function  of  channel 
length  as  illustrated  in  Figure  1g.  It  is  evident  that  the 
nominal  line-defect  density  increases  with  the  channel 
length  and  reaches  a  saturation  level  close  to  15 
20  fiw  (Figure  1g),  which  is  roughly  equal  to  the 
estimated  grain  sizes  for  the  MoS2  films.11  The  results 
presented  in  Figure  1g  serve  as  the  basis  for  our  study 
of  average  GB  effects  in  MoS2  transport  properties. 

For  strategy  2,  we  conducted  transport  measure¬ 
ments  on  different  FET  devices  with  18-^m-long  channel 
length  and  varied  line  defect  densities  and  demon¬ 
strated  that  both  mobility  and  on-state  current  linearly 
decreased  as  a  function  of  the  increase  in-line  defect 
density  from  0.16  to  0.24  /<m_1  (Figure  2a).  The  esti¬ 
mated  mobilities  in  the  presented  range  of  line  defect 
densities  vary  from  4  to  2  cm2/(V  s)  (Figure  2a  inset). 


Field-effect  properties  of  devices  with  varied  channel 
lengths  showed  similar  characteristics  (Figure  2b).  The 
mobility  measurements  show  an  average  decrease  by 
roughly  50%  as  the  channel  length  increases  from  3  to 
~30  /<m,  which  corresponds  to  changes  in  line-defect 
density  from  0.1  to  02  ftm~ The  measurements  also 
show  a  decrease  of  roughly  1  order  ofmagnitude  in  the 
on/off-current  ratios  in  the  same  defect  density  range 
(Figure  2c).  As  the  channel  length  increases,  the  trend 
of  changes  in  the  transport  properties  dosely  follow  the 
changes  in  the  line  defect  density  and  reaches  a  plateau 
region  beyond  20  //m  in  channel  length  where  the  line- 
defect  density  saturates. 

It  is  worth  noting  that  the  study  of  channel  length 
dependence  by  Liu  et  al.  details  the  short  channel 
effects  on  the  transport  properties  of  single  crystal 
MoS2.20,21  They  reveal  that,  for  channel  lengths  larger 
than  ~1  fim,  the  mobility  and  on/off  ratio  are  expected 
to  saturate  at  a  constant  value,  as  the  transport  is 
predominantly  in  the  diffusive  regime.20,21  Given  that 
the  impurity  concentration  and  trapped  charges  on  the 
substrate  are  roughly  uniform  across  the  channel,  it  can 
be  deduced  that  the  observed  changes  in  our  experi¬ 
ments  are  largely  due  to  the  increased  density  of  line- 
defects  introduced  by  the  GBs.  It  is  known  that  grain 
boundaries  play  an  important  role  as  scattering  centers 
that  degrade  the  charge  carrier  mobility  in  materials. 
Previous  calculations  on  the  static  electronic  structures 
of  various  grain  boundaries  in  MoS2  show  that  the 
presence  of  local  states  deep  inside  the  gap  at  Mo 
sites  along  the  grain  boundaries  act  as  charge  traps 
that  undermine  the  electronic  transport  of  MoS2. These 
scattering  centers  reduce  the  mobility  of  charge 
carriers  in  the  samples  in  an  expected  way  while 
the  local  changes  in  the  band  structure  induced  by 
the  defects  degrade  the  gating  process  and  reduce 
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Figure  3.  (a,  b)  Two  representative  devices  for  the  study  of  grain  boundaries'  transport  properties,  (c)  Four  probe  electrical 
transport  transfer  curves  acquired  for  the  devices  presented  in  (a),  (b),  and  pristine  MoS?  samples;  inset  shows  the  logarithmic 
scale  of  the  two  types  of  grain  boundaries  and  pristine  curves,  demonstrating  their  relative  off  current  characteristics.  (d)Two 
probe  electrical  transport  curves  for  tilt  grain  boundaries,  measured  across  U-)  and  parallel  (0 )  to  the  grain  boundaries;  inset 
demonstrates  their  logarithmic  scale,  off  current  characteristics. 


the  on/off-current  ratio.  Additionally,  the  trapped 
substrate  charges  and  local  interaction  variations  among 
the  devices  are  also  of  high  importance,  influencing 
the  measured  transport  properties.  However,  these 
measurements  are  only  approximate  (although  quite 
indicative)  estimations  of  collective  roles  played  by  GBs 
in  large-area  samples.  In  order  to  obtain  more  determi¬ 
nistic  answer  to  the  underlying  mechanisms  explaining 
the  role  of  grain  boundaries  on  electrical  transport 
behaviors,  one  needs  to  explore  such  effects  for  indivi¬ 
dual  GBs. 

Previous  studies  show  that  through  the  regula¬ 
tion  of  the  CVD  parameters,  the  nudeation  in  MoS2 
monolayers  can  be  controlled  to  form  low  density 
crystals  with  different  relative  orientations  (separated 
or  merged  single-crystal  triangles,  Figure  SI  a,  Support¬ 
ing  Information).11  This  allows  for  the  formation  of  a 
variety  of  samples  with  different  GB  types.  Accordingly, 
to  acquire  a  comprehensive  knowledge  of  the  GB 
effects  on  the  electrical  transport  properties  of  MoS2, 
it  is  essential  to  examine  these  GB  variations  individu¬ 
ally.  It  is  known  that  two  main  GB  types,  tilt  and  180° 
twin,  are  present  in  the  aggregates  of  MoS2  single 
crystals.11-14  However,  a  variety  of  these  boundaries 
based  on  differences  in  tilt  angle  and  termination 
elements  at  the  GB  are  possible.  The  previous  mea¬ 
surements  show  a  slightly  higher  conductivity  and 
off-current  for  devices  made  parallel  to  the  180°  twin 


boundaries  but  almost  no  change  in  measurements 
across  the  same  twin  GB.12  They  also  demonstrate  a 
decrease  in  conductivity  across  the  tilt  GBs.  To  better 
understand  the  origin  of  electron  scattering  in  MoS2 
GBs,  we  first  fabricate  multiple  devices  (at  least  four  for 
each  category)  that  contain  only  a  single  GB  (Figure  3a, b) 
and  assess  their  field  effect  characteristic  and  electrical 
transport  properties  and  then  compare  with  devices 
made  on  single  crystals.  In  Figure  3a  (tilt  GB)  and  3b 
(180°  twin  GB),  examples  of  different  MoS2  GBs  exam¬ 
ined  in  our  experiments  are  demonstrated.  This  GB 
classification  presented  here  is  based  on  compari¬ 
sons  with  similarly  merged  triangular  domains  where 
the  GB  structure  was  determined  by  high  resolution 
transmission  microscope  characterization.11-13  It  has 
been  shown  that  the  crystal  orientation  in  each  trian¬ 
gular  domain  is  parallel  to  the  triangles'  sides.11-13 
Therefore,  the  type  of  GB  can  be  estimated  from  the 
domain  mismatch  angle.  However,  in  this  type  of 
analysis,  the  information  about  the  domain  termina¬ 
tion  atoms  cannot  be  extracted.  Nevertheless,  our 
statistical  analysis  of  multiple  devices  studied  in  each 
category  reveals  two  very  distinct  transport  behaviors. 
The  devices  made  on  these  samples  have  a  channel 
width  of  1  // m  and  lengths  ranging  from  2  to  5 //m.  The 
four-probe  conductivity  measurements  for  several  de¬ 
vices  are  shown  in  Figure  3c.  It  is  clear  that  tilt  GBs  are 
the  only  type  with  a  significant  effect  on  the  electronic 
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transport  of  MoS2.  These  results  show  an  average  four- 
probe  field-effect  electron  mobility  of  1 2.1  ±  2.1  cm 2/(V  s) 
and  6  ±  1.1  cm 2/(V  s)  for  pristine  single  crystalline 
samples  and  samples  with  tilt  GBs,  respectively.  This 
demonstrates  a  reduction  in  field-effect  mobility  by  half 
in  the  presence  of  the  tilt  GBs.  In  contrast,  the  devices 
made  on  180°  twin  GBs  show  little  or  no  adverse  effect 
on  electrical  transport  behaviors  of  MoS2.  A  closer  look 
at  the  off-current  and  the  threshold  gating  voltages  for 
these  GBs  also  shows  a  behavior  similar  to  the  pristine 
samples  (Figure  3c  inset),  indicating  negligible  doping 
effects  in  the  materials  as  a  result  of  the  GBs.  Since  tilt 
GBs  seem  to  play  a  more  important  role  in  influencing 
the  transport  property,  we  then  take  a  closer  look  at 
their  behaviors. 

We  perform  two-probe  transport  experiments  to 
investigate  the  electron-transport  behaviors  across 
and  parallel  to  the  tilt  boundaries  in  order  to  delve 
into  the  underlying  physics  of  such  GBs.  The  measure¬ 
ments  for  the  two  representative  devices  as  compared 
to  pristine  samples  are  presented  in  Figure  3d.  Electron 
transport  parallel  to  the  tilt  boundaries  is  more  influ¬ 
enced  by  the  scattering  effects  than  electron  transport 
across  them.  Ourtwo-probe  measurements  on  multiple 
devices  from  each  type  shows  that  the  average  charge 
carrier  field  effect  mobilities  are  42  ±  0.9  cm2/(V  s)  and 
2.0  ±  0.2  cm2/(V  s)  for  transport  across  and  along  tilt  GB, 
respectively,  as  compared  to  5.8  ±  0.3  cm 2/(V  s)  in  the 
pristine  samples.  To  explain  the  observed  behaviors, 
one  could  immediately  refer  to  the  intrinsic  scattering 
properties  of  GBs  related  to  the  changes  in  the  lattice 
potential  as  electrons  move  from  one  grain  to  the  other 
across  the  GB.  However,  this  would  not  explain  the 
behavior  of  those  devices  made  on  parallel  directions 
to  GBs  in  which  two  parallel  conducting  channels  in 
each  grain  could  exist  It  is  also  hard  to  imagine  that  a 
line  of  atomic  defects  parallel  to  a  conductive  channel 
could  dominate  the  charge  carrier  transport  in  a  device 
with  a  1  fim  wide  channel.  This  is  why  we  hypothesize 
that  the  observed  effects  could  be  a  result  of  intrinsic 
and  extrinsic  scattering  centers  accumulated  at  the  tilt 
GBs.  Next,  we  will  provide  further  evidence  to  prove  our 
hypothesis  and  identify  the  sources  of  such  effects. 

The  examination  of  the  off-current  characteristic 
of  the  devices  made  on  varied  orientations  relative  to 
tilt  GBs,  shown  in  the  inset  of  Figure  3d  in  logarithmic 
scale,  indicates  that  the  threshold  gate  voltage  de¬ 
creases  from  pristine  MoS2  to  samples  with  tilt  GBs.  We 
measure  a  gate  threshold  voltage  of  10.3  ±  1.9,  2.5  ± 
2.2,  and  4.0  ±  2.7  V  for  devices  made  on  pristine, 
across,  and  along  the  tiltGBs,  respectively.  This  behavior 
suggests  an  increase  in  n-type  charge  carriers  in  MoS2 
samples  containing  tilt  GBs,  more  so  in  the  devices 
made  parallel  to  them.  One  can  estimate  the  changes 
in  doping  levels  using  n  =  £o£rAVT/te,  where  £o  and  er  are 
the  permittivities  of  free  space  (8.85  x  10~12  F/m)  and 
Si02  (3.9),  respectively;  Al/T  is  the  difference  between 


the  average  threshold  voltage  measured  in  pristine 
and  GB  included  devices,  e  is  the  electron  charge 
1.6  x  10-19C,  and  t  is  the  thickness  of  Si02  layer 
(285  nm).  The  average  increases  in  doping  levels  in 
devices  made  across  and  along  the  tilt  GBs  as  com¬ 
pared  to  pristine  MoS2  are  estimated  to  be  5.9  x  1 011 
and  1.1  x  1012cm-2,  respectively.  Understanding  this 
observed  feature  associated  with  tilt  GBs  in  MoS2  atomic 
layers  could  shed  light  into  their  important  electronic 
degradation  mechanism  for  defect  engineering  in  future 
devices. 

The  origins  of  charge  carrier  types  in  MoS2  are  an 
unresolved  issue,  and  the  domination  of  either  type 
of  charge  carrier  has  been  observed  in  CVD-grown 
samples.22  It  is  known  that  the  sulfur  vacancies  in  this 
material  result  in  deep  midgap  donor  states  below  the 
conduction  band.1323  This  may  be  an  explanation  for 
both  the  nature  of  electronic  transport  degradation 
and  the  doping  levels  seen  in  these  samples.  To  further 
confirm  this,  we  examine  the  temperature-dependent 
transport  behavior  in  these  samples.  The  results  sum¬ 
marized  in  Figure  4a  c  can  be  used  to  identify  the 
nature  of  electronic  states  and  the  role  of  tilt  GBs.  It  is 
evident  from  the  temperature  dependency  of  carrier 
mobilities  that  the  transport  is  dominated  by  hopping 
mechanism  and  mobility  decreases  at  low  tempera¬ 
tures.  We  plot  the  conductivity  as  a  function  of  tem¬ 
perature  using  the  following  relationship 

o  =  o0(T)  exp[  (To/D1/w+1)] 

where  d  is  the  dimensionality  and  equals  2  in  our  case, 
a0(7)  =  AT™,  where  A  is  constant  and  and  m  ~  0.8  1, 
chosen  as  0.8  in  our  analysis.  In  all  cases,  the  transport 
has  a  general  behavior  where  at  low  temperatures, 
below  a  critical  temperature  7*,  the  conductivity  is 
weakly  sensitive  to  temperature  (Figure  4d  f).  How¬ 
ever,  beyond  T*,  the  conductivity  becomes  gradually 
more  responsive  to  temperature  changes.  Additionally, 
the  changes  in  charge  carrier  density  through  elec¬ 
trostatic  gating  show  a  complementary  role  through 
screening  property  changes.  The  observed  behavior 
is  common  in  highly  disordered  systems  and  can  be 
explained  by  Anderson's  model.  This  approach  intro¬ 
duces  random  potential  wells  into  the  periodic  poten¬ 
tial  of  the  material.  The  density  and  depth  of  these 
potential  wells  determine  the  nature  of  disorder,  and 
as  they  are  increased,  the  electronic  wave  function  in 
the  material  becomes  localized.  This  localization  is 
strengthened  as  the  density  and  depth  of  potential 
wells  are  increased.  Free  electrons  are  then  trapped 
in  these  localized  states,  right  below  the  conduction 
band.  At  low  temperatures,  the  electrons  do  not  con¬ 
tribute  to  the  transport  and  play  the  role  of  scattering 
centers.  As  the  temperature  is  increased,  the  electrons 
can  escape  these  states  and  contribute  to  the  current 
in  an  activated  manner.  With  this  picture  in  mind,  one 
can  arrive  to  a  variable  range  and  nearest  neighbor 
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Figure  4.  Temperature  dependent  field  effect  characteristics  measured  using  the  two  probe  configuration  of  (a)  pristine 
samples,  (b)  across  tilt  grain  boundaries,  (c)  along  tilt  grain  boundaries.  Gate  voltage  dependency  of  activation  energy. 
Temperature  dependence  of  conductivity  (<j)  and  variable  range  hopping  (VRH)  at  different  back  gate  voltages,  measured  for 
(d)  pristine,  (e)  across  (1)  the  tilt  grain  boundary,  and  (f)  parallel  (II )  to  the  tilt  grain  boundaries. 


Figure  5.  (a)  Nearest  neighbor  hopping  high  temperature  regime  temperature  coefficient  T0  and  its  gate  voltage  depen 
dence.  (b)  The  variable  range  hopping  low  temperature  regime  temperature  coefficient  T,  and  its  gate  voltage  dependence. 


hopping  transport  mechanism  for  low  and  high  tem¬ 
perature  regimes,  respectively.  The  following  relation¬ 
ships  can  be  therefore  derived  to  express  the 
temperature  dependency  of  the  conductivity  a. 

r<ro«exP[  (r,/nV3] 

f  >  F*  a  «  exp(  To/T) 

The  two  temperature  coefficients  T0  and  T,  are 
important  in  describing  the  nature  of  localized  states. 
Using  these  relationships  one  can  summarize  and 
express  the  observations  in  Figure  4  using  the  gate 
voltage  dependence  of  the  temperature  coefficients 
(Figure  5).  Clearly,  as  the  carrier  densities  are  increased 
electrostatically,  both  temperature  coefficients  de¬ 
crease.  This  is  because  at  high  charge  carrier  densities, 
screening  effects  mitigate  the  role  of  localized  states. 
However,  generally  a  considerable  increase  in  T0  can  be 


seen  as  we  move  from  measurements  on  pristine 
samples  to  samples  with  devices  made  perpendicular 
and  parallel  to  the  tilt  GBs.  The  relationship,  KBT0  = 
explains  the  correlation  between  activation 
energy  and  parameters  of  localization  model,  where  Ke 
is  the  Boltzmann  constant,  is  the  density  of  states  at 
the  Fermi  energy,  and  a  is  the  average  spacing  be¬ 
tween  potential  wells.  Assuming  a  similar  density  of 
states  at  the  Fermi  level,  we  can  attribute  the  changes 
seen  between  measurements  on  the  GB  as  compared 
to  pristine  samples  to  the  decrease  In  average  spacing 
of  potential  wells.  This  decrease  at  Vg  =  20  V  translates 
to  an  increase  in  density  of  potential  wells  by  ~4.9  x 
1012cm-2and  ~13  x  1013cm_2forvertical  and  parallel 
to  tilt  GB  devices  as  compared  to  pristine  samples, 
respectively.  With  this  information,  we  can  now  focus 
on  understanding  the  nature  of  the  scattering  centers 
and  the  random  potential  wells  introduced  by  them. 
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Electron  microscopy  of  the  mechanically  exfoliated 
samples  shows  that  they  typically  contain  a  sufficiently 
large  concentration  of  sulfur  vacancies  that  can  be 
used  to  explain  the  localized  and  activated  transport 
characteristics  in  MoS^23  Several  other  papers,  in  con¬ 
trast,  argue  that  extrinsic  substrate  doping  is  of  higher 
importance  in  terms  of  charge  carrier  density  and 
scattering  properties  in  MoS2.24-27  They  demonstrate 
that  the  nature  of  charge  carriers  in  MoS2  is  closely 
related  to  the  trap  states  in  Si02,  which  can  contribute 
significantly  to  the  conductivity  of  the  material  if  they 
are  close  to  the  conduction  or  valence  band.  These 
observations  provide  two  explanations  for  the  sources 
of  doping  and  scattering  in  the  material.  However,  in 
the  case  of  CVD-grown  MoS2  the  density  of  defects  is 
unsurprisingly  higher  than  samples  prepared  by  ex¬ 
foliating  natural  crystals.13  If  the  source  of  scattering  is 
the  sulfur  vacancies,  then  one  can  only  explain  the 
observed  behavior  by  assuming  a  uniform  distribution 
of  sulfur  vacancies  in  the  CVD  MoS2  and  a  higher 
concentration  accumulated  at  the  GBs.This  is  in  agree¬ 
ment  with  the  previous  observation  of  instability  and 
rapid  decomposition  close  to  MoS2  GBs  at  elevated 
temperatures  (350  K).28  This  explanation  is  also  suitable 
to  explain  the  anomalous  behaviorobserved  in  devices 
parallel  to  the  line  defects,  where,  despite  the  absence 
of  line  defects  in  their  di  rect  path,  they  seem  to  have  an 
even  more  prominent  effect  on  the  transport.  This  way, 
the  devices  with  channels  fabricated  parallel  to  the  line 
defects  contain  higher  average  densities  of  point  defects 
as  compared  to  devices  made  across  GBs  (Figure  S2, 
Supporting  Information).  This  can  be  further  explored 
by  examining  the  T,  gate  voltage  dependence  and  its 
relationship  with  Tq.  We  observe  that,  as  expected,  for  all 
configurations,  T,  decreases  as  Vg  is  increased.  However, 
the  optimum  hopping  distance  results  in  the  following 
relationship  between  7",  and  Tq23 


T\ 


where  £  is  the  localization  length.  The  average  localiza¬ 
tion  length  for  sulfur  vacancies  estimated  using  density 
functional  theory  is  about  6  A.23  Using  the  values 


acquired  in  Figure  5  and  £  for  sulfur  vacancies  we  can 
estimate  the  T,/T 0  ratio  to  test  how  the  model  fits  to  our 
experimental  results.  The  calculated  T,/T0  ratios  acquire 
values  59, 52,  and  45  for  pristine,  across,  and  along  the 
tilt  GBs,  respectively.  The  empirical  values  extracted 
from  Figure  5  for  l/g  =  20  V  are  69,  52, 33,  which  agree 
well  with  the  calculated  values.  As  one  can  see  that  the 
small  discrepancies  between  the  calculated  and  experi¬ 
mental  results  become  significantly  larger  as  the  charge 
carrier  density  is  increased  at  higher  Vg.  This  deviation 
from  the  hopping  model  at  higher  charge  carriers  is 
expected  since  screening  effects  mitigate  the  effects  of 
localized  states  at  higher  electrostatic  gating.  The  agree¬ 
ment  between  the  model  and  empirical  results  at  low 
gate  voltages  signifies  our  hypothesis  for  accumulation 
of  higher  density  of  point  defects  at  the  tilt  GBs. 

CONCLUSION 

In  summary,  a  comprehensive  study  of  the  indivi¬ 
dual  and  collective  role  of  GBs  on  the  transport  proper¬ 
ties  of  MoS2  atomic  layers  was  discussed.  By  examining 
the  collective  and  individual  role  of  line  defects  in 
polycrystalline  MoS2  films,  we  provide  additional  in¬ 
sights  into  the  role  and  mechanism  of  charge  carrier 
scattering.  We  estimate  an  average  decrease  by  half 
in  the  magnitude  of  mobility  in  our  devices  corre¬ 
sponding  to  an  increase  in  the  line  defect  density  by 
0.1  /zm-1.  By  investigating  a  variety  of  GB  types,  we 
observe  that  only  tilt  GBs  have  a  significant  effect  on 
the  mobility  and  transport  characteristics  of  MoS2. 
Further  examination  of  devices  made  along  and  across 
these  tilt  boundaries  reveals  an  anomalous  behavior 
and  a  significant  level  of  n-type  doping.  Low  tempera¬ 
ture  electrical  transport  measurements  confirm  a  hop¬ 
ping  mechanism  for  electron  transport  in  these  samples 
that  can  be  well  explained  with  Anderson  model  for 
disordered  systems.  From  these  observations,  we  con¬ 
clude  that  the  distinct  features  in  tilt  boundaries  are 
more  likely  caused  by  accumulation  of  point  defects  or 
other  structural  defects  at  these  GBs.  This  comprehen¬ 
sive  study  illustrates  the  extent  and  limitations  imposed 
by  GBs  in  the  electronic  transport  properties  of  poly¬ 
crystalline  MoS2  atomic  layers. 


METHODS 

Electron  beam  lithography  was  used  to  fabricate  variable 
channel  length  FET  and  hall  bar  devices.  Photoresist  using  a 
MMA/PMMA  bilayer  resist  process  was  used  for  all  steps.  The 
MoS2  layer  was  patterned  using  a  low  power  ICP  RIE  etch  in  a 
CH4/O2  plasma,  and  source  and  drain  contacts  were  formed  via 
e  beam  evaporated  Ti/Au  (15/85  nm). 

Field  effect  transistor  (FET)  devices  with  variable  channel 
lengths  are  made  by  photolithography  process  using  LOR5B 
as  an  adhesive  layer  and  photoresist  SI  81 3,  mask  aligner  (SUSS 
Mask  Aligner  MJB4),and  02  reactive  ion  etching  (RIE,  Phantom 
III).  The  photoresist  is  removed  by  acetone  and  PG  REMOVER. 
The  electrodes  (Ti/Au,  4  nm/40  nm)  are  deposited  using  an 
e  beam  evaporator.  All  electrical  measurements  are  performed 


in  a  Lakeshore  probe  station  under  vacuum  conditions 
(<10  5  Torr).  These  devices  were  analyzed  using  either  an 
Agilent  B1 500A  or  Keithly  4200  Semiconductor  Device  Analyzer. 

In  all  measurements,  we  estimate  the  charge  carrier  mobility 
in  these  devices  using  the  equation  /a  =  [d/rts/dVbgl[(/./(lVCiVrfs))], 
where  L  and  W  are  the  channel  length  and  width.  The  capaci 
tance  between  the  channel  and  the  back  gate  per  unit  area 
is  estimated  to  be  ~1 2  x  10  4Fm  2  (G  =  e&Jd,  where  e0  =  3.9 
and  d  =  285  nm).The  theoretical  estimate  fbrthe  capacitance  of 
Si02  will  provide  us  witha  lower  limit  estimate  for  the  calculated 
mobilities. 

Micro  Raman  and  photoluminesoence  (PL)  measurements 
were  performed  with  a  WITec  Alpha  300RA  system  using  the 
532  nm  line  of  a  frequency  doubled  Nd:YAG  laser  as  the  excita 
tion  source.  The  spectra  were  measured  in  the  backscattering 
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configuration  using  a  lOOx  objective  and  either  a  600  or 
1800  grooves/mm  grating.  The  spot  size  of  the  laser  was 
—342  nm  resulting  In  an  incident  laser  power  density 
—140  fiVf/pm1.  This  laser  power  was  found  to  be  sufficiently 
low  to  not  to  cause  any  shifting  in  the  both  the  in  plane  and  out 
of  plane  modes  of  the  Raman  signature. 
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